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ABSTRACT 

We study the physical conditions, elemental abundances, and kinematics of the high- velocity clouds 
(HVCs) along the sight lines toward active galaxies HE 0226-4110 and PG 0953+414 using Hubble 
Space Telescope Imaging Spectrograph and Far Ultraviolet Spectroscopic Explorer data. No 21-cm 
H I emission is detected in these clouds, but our observations reveal multiple components of HVC 
absorption in lines of H I, C II, C III, C IV, O VI, Si II, Si III, and Si IV in both directions. We 
investigate whether photoionization by the extragalactic background radiation or by escaping Milky 
Way radiation can explain the observed ionization pattern. We find that photoionization is a good 
explanation for the C II, C III, Si II, and Si III features, but not for the O VI or C IV associated 
with the HVCs, suggesting that two principal phases exist: a warm (T w 10 4 K), photoionized phase 
and a hotter (T = 1 — 3 x 10 5 K), collisionally- ionized phase; the broader line widths of the high ions 
are consistent with this multi-phase hypothesis. The warm HVCs toward HE 0226-4110 have high 
levels of ionization (97-99%), and metallicities ([Z/H] between —0.9 and —0.4) close to those in the 
Magellanic Stream, which lies eleven degrees away on the sky at similar velocities. These HVCs may 
well be stripped fragments of the Stream that have been ionized by the pervading radiation field; they 
have thermal pressures that would place them close to equilibrium in a fully ionized 10 6 K Galactic 
corona with nn = 4 — 9 x 10 -5 cm -3 at 50kpc. The warm HVCs seen at —146 and 125 km s -1 toward 
PG 0953+414 have [Z/H]=— 0.6 ± 0.2 and —0.8 ± 0.2, respectively, suggesting they are not formed 
from purely Galactic material. A mini-survey of the hot, collisionally ionized HVC components seen 
here and in five other sight lines finds that in 11/12 cases, the high ions have kinematics and ionic 
ratios that are consistent with an origin in conductive interfaces, where energy flows into the HVCs 
from a hot surrounding medium and produces O VI and C IV-bearing boundary layers. However, the 
broad absorption wing on the O VI profile toward PG 0953+414 is not completely explained by the 
interface scenario. This feature may be tracing the outflow of hot gas into the Milky Way halo as part 
of a Galactic fountain or wind. 

Subject headings: Galaxy: halo - intergalactic medium - ISM: clouds - ultraviolet: ISM 



1. INTRODUCTION 

Surrounding the Milky Way galaxy lies a popula- 
tion of diffuse gaseous clouds that are not co-rotating 
with the disk. First discovered in 21-cm emission by 
iMuller. Port, fe Raimondl {l963), no single model has 
been able to simultaneously explain the kinematics, 
metallicity, and emission prop erties of these so-called 
high- velocity clouds (HVCs) ijWakker fc van Woerdenl 
I1997|) . The failure to find HVCs distributed throughout 

1 Based on observations from the NASA-CNES-CSA Far Ultra- 
violet Spectroscopic Explorer mission, operated by Johns Hopkins 
University, supported by NASA contract NAS 5-32985, and from 
the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. 

2 E-mail fox@astro.wisc.edu 



other galaxy groups ijPisano et al.1 12004D togethe r with 
HVC detections around M31 l|Braun fc Thilkerl |2004|) 
suggest that HVCs are typically associated with large 
galaxies, rather than the Local Group. Within this 
parad igm, many HVCs (e.g., Complex C. IWakker et alJ 
1999) could represent material leftover from the galax y 
formation process <|OortHl97ft iMaller fc Bullockli2004D . 
or else material stri pped from nearby galaxies (e.g., the 
Magellanic Stream, Gardin er fc Noguchilll99ql . Study- 
ing and understanding HVCs allows us to trace the in- 
teractions and feedback between galaxies and the inter- 
galactic medium. HVCs regulate the rate of star forma- 
tion in our own Galaxy, by providing gaseous material 
from which new stars are born. 

A new phase of gas in HVCs was discovered with 
the detecti ons of high-velocity C IV llSembach et all 
Il995l Il999|) and O VI absorption l|Sembach et all 120001: 
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iMurphv et alJl2000() . A recent survey of extragalactic 
sight lines with the Fa r Ultraviolet Sp ectroscopic Ex- 
plorer (FUSE) satellite i|Sembach et al.ll2003L hereafter 
S03) has found that high-velocity O VI absorption is de- 
tected at >3cr significance 84 times in a survey of 102 
sight lines. Many of these high-velocity O VI features 
are spatially and kinematically associated with known 
21-cm emitting H I structures, for example C omplex C, 
Comp lex A, and the Magellanic Stream (S03, IFox et all 
2004$. There are other cases where O VI HVCs have no 
counterparts detected in H I 21-cm emission. However , 
studies of such highly ioni zed HVCs l|Collins et alJl2f)04 
120051 iGangulv et al.ll2005fl have found absorption in the 
H I Lyman series and in other metal-line tracers of neu- 
tral and weakly ionized gas at the same velocities as the 
high ions. The high ion-low ion connections in HVCs 
are well explained by an arrangement where the hot gas 
arises in the boundary layers between the neu tral HVC 
and the surrounding medium ijFox et al.ll2004fl . 

A separate class of high-velocity O VI absorber exists 
in the form of broad, shallow absorption wings at posi- 
tive velocities. Unlike many of the discrete high veloc- 
ity O VI absorption features seen along complete paths 
through the Galactic halo and corona, the O VI positive- 
velocity wings manifest themselves as continuous ab- 
sorption troughs extending from velocities indicative of 
Galactic halo gas (|^ls,r| < 100 km s -1 ) out to veloci- 
ties of several hundred kms -1 . A majority of positive- 
velocity wings (17/21) are found in the Northern Galactic 
hemisphere, at Z=180-300°and b >30°, even though the 
sample of sight lines is spread fairly evenly around the 
sky (S03). Such a trend may be related to the general 
enhancement of O VI in the Nor thern hemisphere of the 
Galactic halo (|Savage et al.ll2003T) . and to the opening of 
the Local Bubble in this direction l|Welsh et al.1 119991 

In this paper we present high-resolution Space Tele- 
scope Imaging Spectrograph (STIS) observations com- 
plemented with FUSE data to explore the high-velocity 
absorption in one direction showing high-velocity O VI 
components (HE 0226-4110) and one showing a high- 
velocity O VI wing (PG 0953+414). Basic properties of 
the two sight lines under study are found in Table 1. 

The line of sight to Seyfert 1 galaxy HE 0226-4110 
(I = 253.94°, b = -65.78°, z em = 0.495) lies 10.8° from 
the N(H I)= 2.0 x 10 18 cm- 2 (5a) contour of H I 21- 
cm emission from the Magellanic Stream (Figure 1, left 
panel) . The Magellanic Stream is a broad filament of ma- 
terial covering ~1000 square degrees in an arc extending 
from / ss 90°, b « —45°, through t he South Galactic pole 
and around to pos itive latitudes l|Putman et al.l l2003a: 
iBriins et al. ll2005[) . The Stream is believed to be gas 
stripped out of the Magellanic Clou ds by tidal forces as 
the C louds orbit the Milky Way ({Gardiner fe Noguchil 
1996). Six sight lines passing through the Magellanic 
Stream have O VI detections at similar velocities to the 
H I (Figure 15 in S03), suggesting that the Stream con- 
tains (or is surrounded by) a highly ionized component. 

PG 0953+414 is a Seyfert 1 galaxy at z em = 0.239 lying 
at 1= 179.79°, b= +51.71°. Various nearby H I structures 
can be seen in Figure 1 (right panel). This direction 
has the advantage of being very close the anticenter di- 
rection (Z=180°), perpendicular to the direction of LSR 
motion, so the effects of Galactic rotation on absorption 
line profiles are removed. PG 0953+414 is therefore ide- 



ally situated for use as a backlight to study the inflow 
and outflow of gaseous material. 

The structure of this paper is as follows. In §2 we de- 
scribe the observations and data reduction. In §3 our 
spectroscopic measurement and analysis techniques are 
discussed. In §4 and §5 we discuss our observations of the 
high-velocity absorption systems seen toward HE 0226- 
4110 and PG 0953+414, respectively. We discuss pho- 
toionization modeling of the HVCs in §6, and collisional 
ionization modeling of the highly ionized species in §7. 
Section 8 is devoted to understanding the origin of the 
O VI wing absorption seen toward PG 0953+414. The 
principal results of our study are summarized in §9. 

2. OBSERVATIONS AND DATA REDUCTION 

HE 0226-4110 and PG 0953+414 were each observed 
on several occasion s with the E140M eche ll e mode of the 
STIS instrument fWoodgate et all Il998t iKimble et all 
1998), which provides spectra at 7kms _1 resolution 
(FWHM) in the range 1150-1700A, with a few small 
gaps between orders at A >1600A. Details of these ob- 
servations are gi ven in Table 2. Th e data were reduced 
as described in iTripp et"aTl l)2001l) including the two- 
dimensional scattered light correct ion developed by the 
STIS Instrument Definitio n Team ({Landsman fe Bowers! 
1997; Bo wers et al.lll998|) . In gene ral, the wavelength 
calibr ation of STIS data is excellent; iKim Quiiano et all 
(2003) report that the relative wavelength calibration 
is accurate to 0.25-0.5 pixels across the entire spec- 
trum, and the absolute wavelength calibration is accu- 
rate to ~1 pixel (3.2 km s -1 ). The overall velocity uncer- 
tainty in our data was verified by cross-checking the cen- 
ters of Galactic metal-line absorption in O I A1302.169, 
Fe II A1608.451, Si II A1304.370, Si II A1526.707, Ni II 
A1317.217, S II A1253.8 11, and S II A1250 584. 

Th e FUSE satellite i|Moos et al.ll2000t iSahnow et alJ 
l2000j) provided the far-UV spectra of HE 0226-4110 and 
PG 0953+414 (Table 2). For wavelengths below 1000 A, 
we use data from the SiC2A detector segment only. We 
did not combine with data from the SiClB detector seg- 
ment data because this detector has lower sensitivity and 
resolution. For wavelengths between 1000 and 1187 A, 
we use coadded LiFl and LiF2 data, since the resolution 
difference between these two channels is small. The data 
were reduced with the data reduction pipeline CALFUSE 
(v2.1.6 or v2.4.0). Even after the application of the 
pipeline, velocity shifts have to be determined for each 
exposure and each detector segment before coaddition of 
the spe ctra to obtain accur ate wavelength solutions (see 
§3.4 in iWakker et al. 1 120 03). We measured the average 
velocity centroid of all the unsaturated neutral ISM lines 
clearly detected in each FUSE segment (Si II A1020.699, 
Ar I A1048.220, N I A1134.165, 1134.415, 1134.980, and 
numerous H2, Fe II, and O I lines), and then determined 
the offset from the centers of neutral ISM lines seen in 
the STIS data. This "bootstrapping" technique allows 
the kinematics of absorption line profiles observed with 
FUSE and STIS to be reliably compared. All velocities 
quoted in this paper are referenced to the Local Stan- 
dard of Rest (LSR). Note that no velocity difference was 
found between the metal lines and the H2 lines. The off- 
sets for the LiF channels are typically ss lOkrns -1 but 
in the SiC channels they can be as high as ks 40kms . 
We estimate a residual la zero-point error of 5kms _1 



Multi-phase HVCs 



3 



after this process has been applied, since 5kms is the 
dispersion of post-correction velocity centroid measure- 
ments of neutral lines in the FUSE data. This velocity 
error is propagated in calculating the systematic errors 
on our absorption line measurements. 

Continua were fitted to both the FUSE and ST1S data 
in regions approximately 1000 km s -1 wide around each 
absorption line of interest, using regions of the spectra 
judged to be free from absorption. In almost all cases a 
linear continuum could be fitted, since our target active 
galaxies have flat UV spectra generated by non-thermal 
emission processes. 

All FUSE data were rebinned by five pixels to provide 
optimal sampling, for both display and measurement. 
The raw STIS spectra are already optimally sampled, 
so we measure the absorption using the unbinncd data, 
but since the STIS spectra for both HE 0226-41 10 and 
PG 0953+414 are very noisy, we use 2-pixel rebinning for 
display 3 . 

For each of our two studied sight lines we display H I 
21-cm emission spectra, to emphasize the low iV(H I) of 
the highly ionized HVCs under study. For the HE 0226- 
4110 direction these data were obtained and provided by 
R. Morras, using the Villa Elisa telescope, which has 
a beam size of 34'; some residual uncertainty in the 
baseline fit is present in the range ±100 km s^ 1 . For 
PG 09 53+414, we use spe ctra from the Effelsberg tele- 
scope ijWakker et al.ll200l() with a 9'. 7 beam. Both data 
sets have an approximate velocity resolution of lkms" 1 , 
and have been corrected for side lobe contamination. 

3. SPECTROSCOPIC MEASUREMENTS 

We identified the high- velocity component structure in 
each sight line using the high-resolution STIS observa- 
tions of metal-line absorption, particularly in the C II 
A1334.532, Si II A1260.422, Si II A1193.290, and Si III 
A1206.500 lines. We then estimated the velocity range 
covered by each absorption component, and measured 
the equivalent width of the absorption between these 
limits in all spectral lines present in the data. We con- 
verted the measured optical depth in the line as a func- 
tion of velocity into apparent column densities, N a (v), 
using the relation N a (v) = 3.768 x 10 u (/A) _1 T a (u) 
ionscm~ 2 (kms~ 1 ) -1 , where the apparent optical depth 
(AOD) is given by r a (v) = ln[F c (v)/F(v)], / is the os- 
cillat or strength of the transition (taken from iMortonl 
2003), A is the transition wavelength in Angstroms, and 
F(v) and F c (v) are the observed l ine and continuum 
fluxe s at velocity v, respectively (|Savage fc Sembachl 
I1991D . Integrating N a (v) between two chosen velocity 
limits determines the total column density in that veloc- 
ity range, and the optical depth- weighted average values 
of v and y/[v — v) 2 , i.e. the first and second moments 
of T a (v \ give measurements of the line center and line 
width fffemhMh fc Sava g eiri99l. 

The AOD technique requires no prior knowledge of 
the component structure. It will produce good results 
as long as the line profiles are resolved and not satu- 
rated. STIS/E140M, FUSE /HF, and FUSE/SiC can 
resolve lines with b >4, 12, and 15kms _1 , respectively. 

3 Exception: in one case, where we need the highest possible 
resolution to emphasize the component structure, we show the fully 
sampled spectra (§4.3). 



iWakker et all pOO^ show a comparison of the appar- 
ent column density of O VI from both members of the 
doublet toward HE 0226-4110 and PG 0953+414 in the 
velocity range —400 to 400 km s -1 ; this data shows that 
no unresolved saturation affects the zero- velocity (Milky 
Way disk and halo) O VI absorption. By inference, the 
weaker high- velocity O VI absorption, of interest here, is 
also likely to be unsaturated. The O VI lines are likely 
to be resolved since the high temperatures of the regions 
where O VI is found result in lines with thermal widths 
b > 18kms . However, FUSE absorption line profiles 
of the H I Lyman lines may be unresolved, making sat- 
uration hard to detect. As one measures progressively 
weaker Lyman series absorption lines, the /-values of the 
transitions decrease, and eventually the lines will become 
unsaturated. At this point, consecutive Lyman lines will 
give the same value for N a ; this occured for the HVCs 
along our two sight lines with the H I A 920.963, 919.351, 
918.129, and 917.181 lines. To check the column densi- 
ties derived from AOD integrations of these lines, we con- 
struct a Voigt profile model and compare with the data. 
For the low-velocity components, we take the parame- 
ters from the fits to the 21-cm pr ofiles in our two tar- 
get directions (Wa kker et alJl2003|k for the high- velocity 
components, we take the velocities, widths, and column 
densities defined by the AOD integrations. All compo- 
nents are convolved with a Gaussian line spread function 
(with FWHM 25kms _1 for the FUSE /SiC channels). 
Our model also reproduces the D I and O I lines when 
D/H= 2 x 10~ 5 and O/H is solar, assuming they have 
the same width as the H I lines (i.e., non-thermal broad- 
ening dominates). Small adjustments (w 5kms _1 ) are 
needed to the widths of each H I component to produce 
a good fit 4 , but the AOD column densities reproduce the 
depth of absorption well. We thus have confidence in 
the values of iV(H I) derived from the AOD method. A 
full Voigt profile fit of the H I Lyman series absorption 
lines, involving minimizing x 2 f° r & grid of different val- 
ues of log V(H I) and 6(H I), was attempted, but did not 
produce reliable results given the overlapping component 
structure, poor S/N, numerous blends, and degeneracy 
between N and b. 

To check for the possibility of saturation in the strong 
metal lines (C II, C III, Si II, Si III), Voigt profile mod- 
els were constructed in the same manner as for the H I 
lines (using the column densities from the N a (v) inte- 
grations) and compared with the observations. We in- 
vestigated how saturated the observed profiles could be 
by increasing the model component column densities and 
noting whether the quality of the fit decreases. For the 
cases of the C III and Si III lines toward HE 0226-4110, 
saturation is clearly possible in some HVC components, 
resulting in large positive errors (as large as 1 dex) in 
log N. 

4 For PG 0953+414, we found that in addition to the HVCs 
measured here, two extra components were needed to reproduce 
the Lyman series and O I absorption line profiles. A component at 
— 110 kms -1 , with iV(H I)= 1 — 2x 10 17 is needed. This component 
is also seen in the negative-velocity wings of Si II A1260.422 and 
Si III A1206.500, and in H I 21-cm emission at this velocity less 
than 2° away, in HVC Mill (see Figure 1). A second component 
at +55 kms" 1 with N(H I)= 1 - 2 X 10 17 cm~ 2 is needed to fill in 
the H I absorption between the zero-velocity component and the 
positive- velocity HVC. Neither of these components changes the 
H I column density of the HVCs under study in this paper. 



4 



Fox et al. 



The full absorption line measurements are presented 
on a component-by-component basis for the HE 0226- 
4110 system in Table 3, and for the PG 0953+414 sys- 
tem in Table 4. We list the component velocity, width, 
and column density measurements returned by the AOD 
method, the equivalent width, and the observed S/N near 
the line, measured per resolution element in the contin- 
uum. When no absorption is detected in a particular fea- 
ture, we calculate a 3<r upper limit for the ionic column 
density using the error on the equivalent width measure- 
ment, converted to a column density limit assuming a 
linear curve-of-growth 5 . 

Two sets of errors are quoted for each measurement 
of equivalent width and column density. The first is a 
quadrature sum of statistical uncertainty in the count 
rate (Poisson noise) and the continuum placement uncer- 
tainty. The second is the systematic error, found from a 
quadrature addition of fixed pattern noise (f.p.n.) in the 
detectors with the velocity scale errors. We estimated 
the f.p.n. per unbinncd pixel to be 1% for STIS de- 
tectors, 2.5% for FUSE observations of HE 0226-4110, 
and 10% for FUSE observations of PG 0953+414. The 
f.p.n. is lower in the HE 0226-4110 FUSE dataset due 
to the use of focal-plane (FP) splits for some observa- 
tions, which place the spectra in different positions on 
the detector. These errors imply that any feature weaker 
than 1.8 mA (6.3 mA) could be caused by f.p.n. in the 
HE 0226-4110 (PG 0953+414) FUSE data. The veloc- 
ity error (±5kms _1 ) allows for the intrinsic uncertainty 
in the zero point of the post-calibration velocity scale, 
and the uncertainty in the choice of velocities used to 
define a particular high- velocity component. The latter 
uncertainty is important since the different components 
of high-velocity absorption are not always distinct from 
each other and placing the division can be difficult. We 
calculate the effect of this error by changing the velocity 
extrema of each component by ±5kms _1 and observing 
the change in the equivalent width and in N a . 

For each sight line we have prepared three plots show- 
ing the absorption line data (Figures 2, 3, and 4 for 
HE 0226-4110 and Figures 5, 6, and 7 for PG 0953+414). 
First we display (Figs. 2 and 5) the fully-reduced spectra 
for all absorption lines in which high- velocity absorption 
is detected, using gray shading to identify the gas in the 
high-velocity components. We also include H I 21-cm 
emission line spectra of our two sight lines in the upper 
left panels. In Figures 3 and 6 we include stacks of H I 
Lyman series absorption profiles. These figures show the 
presence of high-velocity H I absorption in the Lyman 
series from Ly7 at 972.537 A [n = 1— »4) down to Ly/x 
at 917.181 A (n = 1 — > 13). At wavelengths shortward 
of Ly/x the flux goes to zero. Lya is too broad and Ly/3 
is too polluted by geocoronal emission (O I airglow) to 
show the high-velocity component structure. Geocoronal 
emission in the stronger Lyman lines tends to "fill in" ab- 
sorption in the approximate range —100 to 100 km s" 1 , so 
we display night-only data for the H I lines to minimize 
this effect; at high velocities the geocoronal emission be- 
comes insignificant. Our final method of displaying the 
data is to compare the apparent column density profiles 
of particular lines (Figures 4 and 7) . 

5 The relation used is N = 1.13 X W 17 W X /^ 2 /, where N is 
in cm -2 , W\ is the 3(r equivalent width limit in mA, and A is in A. 



3.1. Blends toward HE 0226-4110 

An IGM absorption line system a t z = 0.34028 is see n 
in Lya, Ly/3, C III, O IV, and O VI (|Lehner et al.ll2005|) . 
Absorption at this redshift by Ne VIII A770.409 would 
blend with high-velocity absorption from the strong line 
of O VI at 1031.926A in the range 150-210kms _1 . We 
used the redshifted absorption data for the weak line in 
the Ne VIII doublet (A780.324) to assess the strength 
of the blend; the weak line shows no absorption at the 
same redshift, with W\(Ne VIII) <9mA (3<r) integrated 
over 60kms _1 (the width of the O VI lines in the IGM 
system). Therefore the strength of the Ne VIII A770.409 
blend must be <18mA (3<r), since /770//78O = 2. The 
total absorption in the range 160-230 km s -1 in the frame 
of the O VI line has W\ = 45 + 5 + 5 mA, more than can 
be accounted for by the blend, so we proceed by treating 
this as a true detection of O VI. 

Absorption from the J = 4 rotational level of molecular 
hydrogen at 1032.156 A blends with high- velocity absorp- 
tion from O VI A1031.926, and absorption from the J = 1 
level at 1038. 156 A blends with high- velocity absorption 
from the O VI A1037.617, if the H2 column densities are 
high enough. We employed a modeling process to de- 
contaminate the data from blending by H2 lines. We 
took all H2 lines in the FUSE data with no major con- 
tamination or continuum placement uncertainties, and 
measured their equivalent widths. We then performed a 
curve-of-growth analysis for each rotational (J) level, and 
derived logA(H 2 ) = 13.72+0.15, 13.99+0.23, <13.89, 
<13.89, and <13.83 in levels J = 0, 1, 2, 3, and 4 respec- 
tively, corresponding to log N (H 2 )«14.34 in total. This 
model is included in the continuum fits for contaminated 
FUSE absorption lines in Figure 2. Examination of the 
O VI A1037.617 panel shows significant contamination of 
the high-velocity absorption by H 2 A1038.156, and so we 
choose to base all measurements of high- velocity O VI on 
the strong O VI line. 

Contamination from redshifted O VI A1037.617 ab- 
sorption associated with the AGN itself (z — 0.49246) 
blends with low- and high- velocity C IV absorption in the 
1548. 195A line. This associated absorber has a multiple- 
component structure and prevents a reliable measure- 
ment of high-velocity C IV A1548.195 absorption. We 
therefore concentrate on the weak line (1550. 770A) for 
measurements of high-velocity C IV absorption in this 
direction. 

3.2. Blends toward PG 0953+4U 

Our measurement of the foreground H 2 absorption in 
the PG 0953+414 direction finds logA(H 2 ) = 13.90, 
14.40, 14.00, 14.12, <14.09 in levels J = 0,1,2,3, and 
4, giving a total log A(H 2 ) = 14.75. As was the case for 
HE 0226-4110, H 2 A1032.156 at 125 km s" 1 in the frame 
of O VI A1031.926 is predicted to be of negligible stength, 
but the J = 1 H 2 line at 1038.156 A (central line depth 
0.27) blends with high- velocity absorption from the weak 
line of O VI at 1037.617A (See Figure 5). We therefore 
adopt measurements of the strong line of O VI only. 

There are few IGM systems ne ar our lines of i nteres t in 
the PG 0953+414 sight line. See lSavage et alJ l|2002fl for 
an investigation of the IGM absorption in this direction. 

4. HIGH- VELOCITY ABSORPTION TOWARD HE 0226-4110 
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Our FUSE and HST/STIS observations of HE 0226- 
4110 reveal high- velocity absorption in lines of H I, C II, 
C III, C IV, O VI, Si II, Si III, and Si IV (Figure 
2). Fe II is detected at low significance. No absorp- 
tion is seen in N I, N II, O I, S II, S III, P II, Ar I, 
and Fe III. Multi-component structure is clearly seen 
in the HVCs in the Si II, Si III, C II, and C IV pro- 
files observed at 7 kms -1 resolution. Based on compo- 
nent fits to these lines, we identify four positive-velocity 
components centered at 99, 148, 175, and 193kms~ 1 , 
covering the ranges 80-125, 125-160, 160-185, and 185- 
230 km s , and shown in Figure 4 with the colors pur- 
ple, blue, green, and yellow respectively. Absorption in 
the 99 kms" 1 component is distinct from the Galactic 
absorption centered at — 10 kms -1 , so we treat this com- 
ponent as an HVC, even though the blueshifted parts of 
this component lie outside the formal definition of high- 
velocity (\vlsr\ > 100 kms" 1 ). 

Inspection of Figure 4 reveals the HVC component 
structure to be consistent among lines of different ion- 
ization states. We show the STIS data at full resolution 
here so that no kinematic information is lost. Tracers of 
neutral (C II, Si II), weakly ionized (Si III), and highly 
ionized gas (C IV, Si IV, O VI) show absorption over 
the same velocity range, and with a similar component 
structure. Among the high ions, the 99 (purple), 175 
(green), and 193 kms" 1 (yellow) components are clearest 
in O VI, C IV, and Si IV, respectively. This indicates that 
the high-ion/high-ion ratios differ between components. 
There are further differences between components in the 
low-ion/high-ion ratios; for example, N(Si II)/iV(Si IV) 
changes from > 2.0 at 145 kms -1 to 1.0 at 193 kms" 1 . 

The unblended O VI absorption in the 99 kms -1 com- 
ponent is broader than the corresponding lower ioniza- 
tion component; we measure 6(0 VI)=33±2kms -1 and 
b{C II)=16±3kms -1 . The fact that high-velocity ab- 
sorption is seen over a similar velocity extent in both neu- 
tral, weakly ionized, and highly ionized species indicates 
that common structures contain all these ions. However, 
the line widths imply that the neutral and highly ion- 
ized species are not co-spatial, and so multiple phases 
are required in each component. 

The H I Lyman lines (Figure 3) show strong high- 
velocity absorption extending out to 230 kms -1 , the 
same extreme velocity as is seen in the O VI profile. We 
note that the H I column densities in the four HVCs 
toward HE 0226-4110 all lie between 10 16 ' 2 and 10 16 7 . 
Since H I clouds become optically thick to hydrogen ion- 
izing radiation (A < 912 A) at logiV(H I) « 17.2, these 
HVCs are optically thin. Consequently, there should be 
little shielding of the HVC core, and we expect the ion- 
izing radiation field to be relatively constant within the 
clouds. 

5. HIGH- VELOCITY ABSORPTION TOWARD PG 0953+414 

FUSE and STIS data of PG 0953+414 reveal high- 
velocity absorption in H I, C II, C III, C IV, Si III, Si IV, 
and O VI in the range 80-200 kms" 1 , together with a 
narrow —146 kms" 1 HVC seen in the lines of C II, Si II, 
and Si III (Figure 5, Table 4). 

In the positive- velocity HVC, a single, narrow (b = 
10 kms" 1 ) neutral component at 124 kms" 1 is clearly 
seen in the C II profile (Figure 7, blue shading), and at 
moderate significance (3. Oct) in the O I profile, though 



no coincident Si II absorption is detected at significant 
levels. The C III and Si III profiles (green), tracing gas 
of intermediate ionization level, also show a component 
centered near 125 kms" 1 , though the absorption in these 
species is much broader (b — 40 kms" 1 ); the C III profile 
also contains a wing extending out to 200 kms -1 . 

The highly ionized species C IV and Si IV (red) show 
weak absorption components centered near 150 kms" 1 at 
detection significances of 2.7 and 2.3cr, respectively, dis- 
placed in velocity from the C III centroid by «25kms" 1 . 
The O VI profile in this direction shows a distinctive 
positive- velocity wing extending from 75 kms -1 out to 
200 kms" 1 . The wing is qualitatively different from 
the absorption profiles in neutral species, weakly ionized 
species, and the highly ionized species C IV and Si IV. 
Interestingly, there is more O VI absorption at velocities 
where H I is strongest (80-125 kms" 1 ) than where C IV 
is strongest (135-200 kms" 1 ); the N(C IV)/iV(0 VI) ra- 
tio therefore changes with velocity over the range of the 
high-velocity absorption. 

The negative- velocity HVC at —146 kms" 1 is seen in 
lines of C II, Si II, Si III, and Si IV, all of which show 
b- values in the range 12-15 kms -1 . Although the H I ab- 
sorption in this component appears to be blended with 
low-velocity absorption, we conclude high- velocity H I is 
present since our Voigt model without the —146 kms -1 
component clearly does not fit the data (Figure 6). The 
— 146 kms" 1 HVC has the interesting property of no as- 
sociated O VI or C IV absorption. 

6. PHOTOIONIZATION MODELLING: THE CLOUD CORES 

We used the 1-dimensional plane-parallel p hotoion- 
ization code CLOUDY (v94.00; iFerland et al7ll998h to 
model the HVCs as uniform density slabs exposed to an 
ionizing radiation field. The purpose of this exercise is 
to assess whether photoionization is a viable ionization 
mechanism, and if so, to provide information on the phys- 
ical conditions and elemental abundances in the clouds. 
We consider two incident radiation fields: the extragalac- 
tic background (EGB) and the radiation field escaping 
from the Milky Way (MW). 

6.1. Extragalactic Background Radiation 

The extragalactic background (EGB) radiation con- 
tains enough ionizing photons to mai ntain the lev- 
els of i onization se e n in m any HVCs i Scmbach et_aD 
2m2 iTripp et all 120031 ICollins et alJ 12004 120051 
Ganguly et alJl2005fl. The E GB field we use is taken 
from lHaardt fe Madaul l)1996|) . who calculate the mean 
intensity J„ as a function of wavelength at various red- 
shifts by following the radiative transfer of a model QSO 
UV spectrum through the IGM (see their Figure 5). The 
EGB field at z = has a mean intensity at the Lyman 
Limit (912 A) J v = 10" 23 ergcm" 2 s -1 Hz -1 sr -1 , and is 
shown as the blue line on Figure 8. 

A CLOUDY calculation requires a geometry in which 
the radiation field is plane-parallel and incident from one 
side of the cloud. For the EGB, the true field is isotropic, 
but since the HVCs under study are optically thin to Ly- 
man continuum radiation, the direction of incident ion- 
izing photons is unimportant. We therefore collapse the 
EGB field into a one-dimensional flux with the same en- 
ergy density as the isotropic field (F v — AttJ^). When 
one divides F v by the photon energy hv and integrates 
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from to 912 A, the EGB field at z = is found to 
contain 10 hydrogen-ionizing (hard) photons cm~ 2 s _1 . 

6.2. Escaping Milky Way Radiation 

To determine the flux of escaping Galactic radia- 
tion, we first estimate the ionizing spectrum in the 
disk, and then consider what fraction of photons es- 
cape into the halo. The basic th eory is outlined in 
IBland-Hawthorn fc Malonevl (|1999|) although we now 
present an extended model. First, the 90-912 A region is 
dominate d by O-B stars and these are confin ed to spiral 
arms (see Bland-Hawthorn & Maloncy 2002) . The unex- 
tinguished O-B star radiation field is normalized to the 
distribution of O stars ne a r the solar circle compiled by 
iVacca. Garmanv. fc Shulil (fj.996). The mean opacity in 
the 90-912 A window is taken to be r=2.8, implying 6% 
of the hard O-B star photons escape into the halo nor- 
mal to the disk. This escape fraction is established from 
the intensity of Ha emission from HVCs with known dis- 
tance bounds llBland-Hawthorn fc Malonevl 119991 120021: 
iPutman et al.l2003b|) . assuming that the Ha is photoion- 
ized. 

Secondly, the disk spectrum in the range 912-2460A is 
derived fro m the observed spectrum in the s olar neigh- 
bourhood ijMezeer. Mathis. k, PanagiaH T982). Here we 
make a key assumption that the observed field at the 
Sun's position has already been extinguished with a mean 
912-2460A opacity of r=1.5. This takes account of the 
fact that some of the field will have originated close to 
the O-B stars, and some from stars which have moved 
away from the main star forming regions. We adopt the 
same opacity normal to the disk such that 22% of the 
soft UV photons escape into the halo. 

The unextinguished 90-2460 A ionizing spectrum in 
the disk is shown as the uppermost green line in Fig- 
ure 8. The extinguished spectra at distances of 10, 50, 
and 100 kpc along the PG 0953+414 sight line are then 
shown with the lower curves. The geometry of both the 
hard and soft MW radiation fields is illustrated in Fig. 
9, where the contours are lines of constant ionizing flux 
$, and the numbers given represent the logarithm of $ 
in units of photons cm~ 2 s _1 . The Galactic hard radia- 
tion field dominates the EGB hard radiation inside the 
log$ = 4 contour, which occurs at approximately 40 kpc 
for directions \b\ < 30°, 100 kpc for 30 < |b| < 60°, and 
180kpc for directions \b\ > 60°. We do not including 
ionizing photons from hot white dwarfs, planetary neb- 
ulae, and Galactic soft X-ray sources, which may affec t 
the MW radiation field ijBregman fc Harringtonl fi986L 
or ionizing photons from the Magellanic Clouds. 

6.3. General Properties of Photoionization models: 
Galactic vs Extragalactic 

In Figure 10 we plot N{C m)/JV(C II), 
N(C IV)/N(C II), N(Si III) /JV(Si II), and 
Y(Si IV)/Y(Si II) against U (= n 7 /n H ) for three 
photoionization cases: pure MW, pure EGB, and 
MW+EGB at 50 kpc toward HE 0226-4110. The 
differences between these curves are caused by the EGB 
being harder than the MW field. At logU < -2.0, this 
leads to higher iV(Si IV)/iV(Si II) and iV(C W)/N(C II) 
ratios in the EGB. At logZJ > -2, the EGB begins 
to ionize silicon and carbon beyond Si IV and C IV, 



resulting in the turnover in the EGB curves shown 
in Figure 10; the MW field produces Si IV and C IV 
as the dominant ionization stages at these ionization 
parameters. The results are robust to increasing the 
number of non-hydrogen-ionizing photons. 

6.4. Modelling the HVCs toward HE 0226-4110 and 
PG 0953+4 U 

We explicitly modeled each HVC absorption compo- 
nent toward HE 0226-4110 and PG 0953+414 by com- 
paring the CLOUDY model predictions with the ob- 
served ionic column densities of carbon and silicon (sum- 
marized in Table 5), and determining the best-fit values 
of logCZ and metallicity [Z/R] 6 . We assume all the H I 
exists in the same phase as the low ions. This is re- 
peated for assumed HVC distances of 10, 50, and 100 kpc, 
with the Galactic radiation field scaled appropriately us- 
ing Figure 9. At each distance we add the uniform EGB 
contribution and enter the combined field into CLOUDY 
(we do not account for the attenuation of EGB photons 
passing through the H I disk of the Milky Way). Table 6 
shows the flux of the radiation field with distance along 
the sight lines to HE 0226-4110 and PG 0953+414. 

6.4.1. Low versus High Ions 

For all HVCs toward both HE 0226-4110 and 
PG 0953+414, and for both EGB and MW radiation 
fields, we found we could reproduce the column densities 
of the singly- and doubly-ionized species (i.e., C II, C III, 
Si II, Si III), but not the triply- ionized species (C IV and 
Si IV) or the O VI. In Figure 11 we show the results from 
an example run, (HVC2 toward PG 0953+414), for two 
different radiation field cases (EGB and MW at 10 kpc). 
The observations are shown as the larger data points, 
with la error bars. For either MW or EGB fiels, the 
observed levels of Si IV, C IV, and O VI are significantly 
(at least 1 dex) higher than the photoionization model 
predicts at the value of U which explains the singly and 
doubly ionized species. The underproduction of the high 
ions is worse for the MW case. We attempted separate 
high-ion only models, but found that we could not simul- 
taneously reproduce the observed iV(Si IV), N(C IV), 
and 7V(0 VI) with one value of U. These results strongly 
suggest that Si IV, C IV, and certainly O VI are not pho- 
toionized in these HVCs. The difficulty in photoionizing 
O VI has been reported elsewhere, both in the modelling 
of observed HVC ionization pat terns (S03, Collin s~et alJ 
I2004L 1200.4 iGangulv et alJl2005() and in th eoretical work 
on the structure of phqtoion ized clouds l)Kepner et alJ 
1999| IGnat fc Sternbe^l200l . 

Given that the species which appear most consistent 
with photoionization are C II, C III, Si II, and Si III, we 
derive our best- fit values of log[/ and metallicity [Z /H] 
(assuming a solar Si/C ratio) for each component by min- 
imizing the x 2 statistic using data from these four ions 
only. By investigating the depth of the \ 2 minima, we 
find typical errors on logf/ and [Z /~H] of ±0.2 dex (95% 

6 [Z /U]=log (N Z /N H ) - log(N z /N H ) e . The relevant 
solar abundances we adopt throughout this section, us- 
ing the definition Ax=n(X) /n(H), are: log Ac = —3.61 
lAllende Prieto Lambert, fc Asplundl 120021): log Am = -4.07 
IHoTwegeitooOlfl: log A r. = -3.34 lAsplund et alJI2004fl: logAg; = 
-4.46 IHolw cgcr 2001); logA Pe = -4.55 IHolwegeij|200ll) . 
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c.L). By combining [Z /H] and \ogU with the tempera- 
ture and ionization breakdown for each ion returned by 
the CLOUDY model, we derive the following quantities 
for the HVCs (see Table 7): gas density, physical size, 
thermal pressure, hydrogen ionization fraction, and pre- 
dicted intensity Ha emission. 

For a given component, the closer the clouds are as- 
sumed to be to the Milky Way, the higher their density, 
pressure, and predicted Ha intensity, and the smaller 
their physical depth. The temperature of the photoion- 
ized clouds remains confined to a range between 1.0- 
2.2xl0 4 K. 

6.4.2. HVCs toward HE 0226-4110 

For the four HVC components toward HE 0226-4110, 
we find overall metallicities of —0.6 ± 0.2, —0.4 ± 0.2, 
—0.9 ±0.3, and —0.6 ±0.2; these derived abundances are 
essentially independent of the assumed distance to the 
clouds. Elemental abundance studies of the Magellanic 
Stream, lying eleven degrees away from the HE 0226- 
4110 sigh t line, have found [Si/H l=-0.8±0.2 in the Lead- 
ing Arm l|Sembach et al.ll2001bft. and fSi II/H1> -1.5 in 
the main body of the stream IjCibson e"t~a l. 2000). Mag- 
ellanic Stream abundances are similar to those in the 
Small Magellanic Cloud 7 , ch aracterized by [Si /H] ,9 m o = 
-0.51, [C/B]smc = -0.66; lIRussell fe Pop^IT^ cor- 
rected to our updated solar abundance scheme). Since 
our derived values for the silicon abundances in the 
HVCs toward HE 0226-4110 are suggestive of a Mag- 
ellanic Stream origin, we are able to place the HVCs at 
«50kpc. Looking at the inferred densities, pressures, 
and cloud sizes of the MW±EGB 50 kpc models of the 
HE 0226-4110 clouds in Table 7, we find that the four 
HVCs have derived densities in the range logriH=— 2.3 
to —2.1, physical sizes of 39 to 120 pc, thermal pressures 
of 95-190 cm~ 3 K 8 , and hydrogen ionization fractions of 
98%. If we assume spherical symmetry for the clouds, 
they have masses in the range 10-250 M Q and angular 
sizes on the sky of several arcminutes. Other geome- 
tries could result in substantially larger values for both 
the mass and angular size. One explanation is that 
the clouds trace material tidally stripped from the main 
(21-cm) body of the Magellanic Stream, which are then 
ionized by the pervading MW radiation field. 

If we treat the detection of Fe II in HVC2 at 2.4a signif- 
icance as real, then the CLOUDY model returns [Fe/H]« 
+0.6 in this component, whe reas carbon and silicon give 
[Z/H] = -0.4. We note that ILu et all lfl998) concluded 
that dust was present in the Magellanic Stream, since 
the measured S/Fe ratio of eight times the solar value 
in HVC 287.5±22.5±240 requires much of the iron to 
be depleted out of the gas phase. Unless the elemen- 
tal abundances are markedly different in the core toward 
HVC 287.5+22.5+240 than in the HVCs under study 
here, we conclude that the iron feature in HVC2 is likely 
n ot real; it cou l d be d ue to fixed pattern noise. 

iCollins et alJ l)2005f ) report a size of 5.9 kpc for the 
same HVC system toward HE 0226-4110 that we model 

7 In the Gardiner & Noeuchi ( 1996) tidal model, all the gas in 
the Magellanic Stream comes out of the SMC. 

8 The main body of the Stream cont ains clumps with higher ther- 
mal pressure: ISembach etail {2001b) measure P/k ~ 360 cm 3 K 
in the leading-arm. 



here, based on different photoionization models. The dif- 
ference between their result and ours (HVC sizes of 39- 
120 pc if the clouds are at 50 kpc) is partly due to their 
modeling the system as a single component, whereas 
we use four components, with smaller individual column 
densities. A second reason is metallicity: we derive val- 
ues hi gher than the 0.1 solar assumed by Collins et al. 
( 2005) , resulting in our finding lower ionization parame- 
ters, higher densities, and smaller cloud sizes. 

6.4.3. HVCs toward PC 0953+4 U 

The metallicity of the 125 km s" 1 HVC toward 
PG 0953+414 is [Z/H] = -0.8±0.2 independent of dis- 
tance, suggesting it is not a fountain cloud, which would 
be expected to have solar metallicity or above, but rather 
a cloud of non-Galactic origin that happens to be moving 
away from the LSR. Our modeling finds the cloud to be 
99-100% ionized. The detection of O I at apparent 3. Oct 
significance in HVC2 implies [O I/H I]«[O/H]=+0.6±0.2 
if the feature is real. Since the carbon and silicon lines 
suggest metallicities lower by «1.6dex, this O I detec- 
tion is likely spurious. We have no information on the 
distance to this absorber. 

The -146 km s" 1 HVC has a metallicity of -0.6±0.2. 
The Si II /Si III and Si II/C II column density ratios (and 
overall line strengths) in the — 150kms _1 HVC toward 
PG 0953+414 are similar to those of the +100 km s" 1 
H VC in the spectrum of P G 1116+215 recently presented 
by iGangulv et alJ l)2005|) . The properties of both of 
these HVCs, including these low-ion ratios, the absence 
of C IV and O VI, and the approximate H I columns, are 
highly reminiscent of the Virgo Cluster L yq cloud in the 
spectr um of 3C 273.0 at z a bs = 0.00530 ijSembach et alJ 
2001a; Tri pp et alJ I2002D . Tripp et al. find that the 
3C 273.0 Lya cloud is quite small, ~70pc, of the same 
order as the sizes w e derive for the — 146 k m s" 1 HVC 
(80 pc if at 100 kpc). iTumlinson et alJ l|2005[) report two 
Lya absorbers toward QSO PG 1211+143, each lying 
within 150 kpc of a potential host galaxy, which also have 
ionization properties similar to those of Galactic highly 
ionized HVCs. 

6.5. HVC Thermal Pressure 

We compare our derived HVC pressures with the es- 
timated pressure of an isothermal Galactic corona, as 
a function of distance. The coronal pressure has been 
calculated bv iSternberg. McKee. fe Wolfird (|2002) for a 
2 x 10 6 K hydrostatic halo using a two-component (disk 
and dark-matter) Galactic potential. Sternberg et al. 
find P H i]\i/k «[1000, 250, 100] cm" 3 K at r =[10, 50, 
100] kpc. By averaging the results of the six HVCs we 
model, we find a typical P/k of [530, 140, 50] cm~ 3 K 
if the clouds are at these same distances. Ignoring the 
difference between the Sun-HVC and Galactic Center- 
HVC distances, our clouds would therefore be slightly 
underpressurized in a 2 x 10 6 K corona, but would be 
close to thermal pressure equilibrium in a 1 x 10 6 K 
corona. Such a medium is suspected f rom UV (S03, 
Fox et al.ll200l ICollins et a,l J EoM 120051 IGangulv et all 
2005) and X-ray high-ion absorption ijFang et al.ll2003t 
McK ernan. Yaaoob. fc Revnoldsl2004j) studies, and from 
pulsar dispersion measures, X-ray emission, and proper- 
ties of the Magellanic Stream (see iQuilis fe Moot3I2001L 
and references therein). At 50 kpc, the density of a fully 
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ionized 10 6 K corona would be 7Jh = 4 - 



9 x 10~ 5 cm" 3 



if 



thermal pressure equilibrium holds in the HVCs toward 
HE 0226-4110, assumed to be part of the Magellanic 
Stream. Magnetic fields and cosmic rays may also pro- 
vide pressure confinement, but the similarity of thermal 
pressures between a hot corona and these HVCs is strik- 
ing nonetheless. 

7. COLLISIONAL IONIZATION: THE CLOUD BOUNDARIES 

Since photoionization cannot explain the high ions, 
we now investigate collisional ionization models, 
specifically: collisional io nization equilibrium (CIE; 
iSutherland fc pop italll99 ffi. magnet ized conductive in- 
terfaces fCIs: iBgrkgw^kletjjjl lOOOl) . turbulent mixing 
layers fT MLs:|S hivhiet^L|[T993 j). radiativelv cooling hot 
gas (R,C: | Ecj gar fc Chevalierll l986t) . and shock ionization 
(SI; IDonita fc Sutherland! 11 A brief description of 

each of these models, and t he appropr i ate ra nge of their 
parameter space, is given in lFox et al.l ((2004) . 

To increase the sample of observations over the high- 
ion detections in HVCs toward HE 0226-4110 and 
PG 0953+414, we produce a compilation of all mea- 
surements of highly ionized HVCs where Si IV, C IV, 
and N V data exist iIFox et alJl200l iCollins et al1l2004t 
iGaneulv et al.ll2005|) . This compilation is presented in 
Table 8. In cases where only a limit could be measured 
for AT (Si IV), N(C IV) or V(N V), we present a corre- 
sponding 3cr upper limit to the ratio of that ion to O VI. 
We calculate the mean and standard deviation of the high 
ion ratios in these fourteen HVCs, and compare these to 
the integrated aver ages measured in the Galactic halo 
ijZsargo et al.ll2003|) . 

There is an enormous range in N(H Tj/N(0 VI) ob- 
served in HVCs, from ~ 10 6 in Complex C to ~ 10 2 
in components toward PG 0953+414 and PKS 2155- 
304. The variation in the ratio is primarily due to vari- 
ations in V(H I), whereas V(0 VI) is always found 
to be within ldex of 13.0. This represents a cru- 
cial piece of evidence in understanding the structure 
of HVCs: no matter how much neutral gas there is, 
some process regulates the production of highly ionized 
gas and always produces a similar amount of O VI. A 
corresponding dispersion in the N(H I)/V(0 VI) ratio 
in low redshift intergalactic O VI absorbers has been 
noted <|Tripp. Savage, fc .Ienkinsll2000l iDanforth fc ShuTll 
2005J. 

Our observations are compared with the model 
predictions in Figure 12, which contains plots of 
AT(C IV)/V(0 VI) vs V(N V)/N(0 VI) (left panel), and 
V(Si IV)/V(0 VI) vs AT(N V)/V(0 VI) (right panel). 
All measurements of high-ion column density ratios in 
HVC components (given in Table 8) are shown in red, up- 
per limits are shown in purple, and the Gal actic halo av- 
erage is shown in green (|Zsarg6 et al.ll2003l . The shaded 
regions represent the range of parameter space allowed 
by the various models. The blue arrows show the effect 
of correcting the model predictions to a scheme of SMC 
abundances l)R,ussell fc Dopitall992|L assuming tha t oxy- 
gen is the dominant coolant (see IFox et al.l 12004. for a 
description of the correction process). A summary of the 
ability of each model to reproduce the observations is 
summarized in Table 9. 

The shock ionization model, when calculated using a 
shock speed of 200 km s _1 (closest to the LSR velocities of 



the HVCs toward HE 0226-4110 and PG 0953+414) sys- 
tematically underestimates N(C IV)/V(0 VI). A much 
higher velocity shock (~ 500 kins ) is needed to ex- 
plain the N(C IV)/V(0 VI) ratio by shock ionization, 
which makes this scenario unlikely, though we cannot 
rule it out, since a HVC falling in at 200 km s -1 against 
a 300 km s^ 1 Galactic outflow would produce a relative 
velocity of ~ 500 km s -1 , and hence a higher Mach num- 
ber. The region of ratio-ratio space in Figure 12 covered 
by the shock model prediction allows for the possibility 
of shock velocities of up to 500 km s -1 . 

Conductive interfaces can explain the 
N(C W)/N(0 VI) and V(N V)/V(0 VI) ratios in al- 
most all cases, but cannot explain the N(Si 1V)/N(0 VI) 
observations. However, less weight should be given to 
the ratios involving V( Si IV). Ionizing photons produced 
in th e interface itself (Scmbach & Savage Hl992t ISlavinl 
120041) might inflate N(Si IV), rendering the purely 
collisional ionization model predictions invalid. In par- 
ticular, the interface should be a strong emitter of He II 
A304 radiation, capable o f ionizing Si IV but not C IV 
llSembach fc Sava"gel 119921: iSavaee. Sembach. fc CardeTIil 
119941) . Such ionizing photons are not accounted for in 
our CLOUDY photoionization models of §6. In the 
case where all ratios involving V(Si IV) are ignored, 
the conductive interface model is the most successful 
at explaining the remaining ratios (C IV/O VI and 
N V/O VI) in HVCs (and the average Galactic halo 
value). The interface model is also consistent with 
the observed HVC absorption kinematics, since the 
observed high ion-low ion offsets are of the same order 
(0-25kms _1 ) as would be expe cted in an evaporative 
flow off a conductive interface (Bohring er fc Hartauistl 

To further our understanding of conductive inter- 
faces, the effects of HVC motion on the conduction 
across the interface need to be explored theoretically. 
Numerical simulations of the interact ion between in- 
falling H VCs and surroun ding media ijQuilis fc Moord 
120011 iMurrav fc Linl 120041) indicate that cometary fea- 
tures, turbulence, and bow shocks are all expected to 
occur around HVCs. More advanced models of con- 
ductive interfaces would also allow for a full treatment 
of the ionizing radiation produced within the inter- 
face i tself (the diffuse field: iFerrara fc Fieldlll99l ISlavinl 
|2004[) . Such a treatment would be able to make pre- 
dictions on the relationship betwe en HVC O VI ab- 
sorption and HVC Ha emission dWei ner fc Williams 
[T9MlTufte. Reynolds, fc Haffnerll998HTufte et al.l2n02t 
iPutman et al.ll2003blK If a significant fraction of the Ha 
emission from HVCs arises from conductive interfaces, 
rather than by photoionization from escaping Galactic 
radiation, then the true escape fraction of Galactic ion- 
izing photons may be less than the 6% (normal to disk) 
discussed in §6. Further Ha studies of HVCs will help 
elucidate this issue. 

8. POSITIVE- VELOCITY WINGS 

The O VI profile toward PG 0953+414, with its smooth 
positive- velocity wing rather than discrete component, is 
qualitatively different from the C II, C III, C IV, Si III, 
and Si IV absorption line profiles in this direction. The 
C IV and Si IV HVCs, although observed at low sig- 
nificance, each appear to show a broad component cen- 
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tered near 150kms _1 . One interpretation of the data is 
that the 125 km s -1 HVC (seen in C II) is surrounded 
by highly ionized interface layers, which contain the ions 
C IV and Si IV. This explanation is consistent with the 
C IV/Si IV ratio, the C IV-C II offset, and our explana- 
tion for the high ions in the HE 0226-4110 clouds and 
other HVCs, but poses the question: where is the O VI in 
this interface? We show using the blue line in the upper- 
right panel of Figure 7 that the interface O VI could well 
be hiding underneath the wing. This line shows the con- 
tribution to the O VI profile of an interface component 
having the same width as the C IV HVC and a typi- 
cal N(C TV)/N(0 VI) interface ratio (of 1). As can be 
seen, such a component can explain all the O VI absorp- 
tion between 140 and 200 km s -1 , but not the remaining 
high-velocity O VI between 75 and l^Okms^ 1 . 

The O VI seen at 75-140 km s" 1 toward PC 0953+414 
has no accompanying absorption seen in C IV or Si IV: we 
find logiV(0 VI)=13.67±0.08, log AT (C IV)<13.22, and 
log N(S\ IV)<12.49 over the same velocity range. Assum- 
ing solar C/O and Si/O ratios and CIE conditions apply, 
the C IV/O VI and Si IV/O VI ratios imply logT > 5.31 
and > 5.26, respectivel y. This hot gas may be tracing a 
hot Galactic Fountain ^Shapiro fe Field! Il976l iBree man 
IT98r3lNorman fc Ikeuchll989UHouck fc Bregmanll990|) : 
a circulation of interstellar material in which supernovae- 
heated gas rises buoyantly into the halo, cools, and then 
rains back toward the disk. In a hot fountain, little O VI 
is expected in the initial 10 6 K flow, but as the gas cools, 
higher stages of ionization will recombine into +5 as the 
flow drives outward, even though gas at 3 x 10 5 K (where 
O VI is formed by collisions) is not hot enough to drive 
the flow. The model is attractive for explaining positive- 
velocity wings since it naturally explains why the wings 
are always seen to be redshifted; no negative- velocity (in- 
falling) wings are seen in O VI profiles. 

The Galactic Fountain model has had some success ex- 
plaining the kinematics and metallicites of the interme- 
diate velocity clouds (IVCs), defin ed to be at \vlsr\ = 
30 - 90 km s" 1 (jRichter et al.l200H) . One could speculate 
that hot outflowing fountain gas is detected as positive- 
velocity wings, and cooler inflowing gas is detected in 
the form of neutral IVCs. In this case, the outflowing 
gas must be dense enough to cool in the time it takes to 
rise out into the halo, otherwise the H I IVCs will never 
condense. Since the instantaneous ra diative cooling time 
is t co o i = 2.4 x 10 8 nlJr 6 Al^ 9 yr l|Houck fc Bregmanl 
1990), where A_ 2 i 9 is the cooling function in units of its 
value at 10 6 K (1.26 x 10~ 22 erg cm 3 s" 1 ), T 6 = T/10 6 K, 
and n_3 = n/10~ 3 cm~ 3 , we find that a 10 6 K fountain 
flow moving at 100 km s -1 must have a density n_3 w 9 
to radiate half its internal energy in a sound crossing 
time (28Myr for a scale height of 2.5 kpc, appropriate 
for an isothermal 10 6 K halo), which approximates the 
time taken for the gas to fill the halo. 

A related possibility is that the wings trace an es- 
caping Galactic wind, energetic enough to escape the 
Galactic potential well. Cosmic rays could play a role 
in dri ving such a Galactic wind and support ing high- 
z gas (|Breitschwerdt. McKenzie. fc Volklll99H) . O VI- 
containing outflows are pre dicted by hydrodynamical 
simulations of galaxy wind s ijMac Low fc Ferraralll999l 
iStrickland fc Stevens! I2000|) . and have been detected in 



starburst galaxies NGC 1705 llHeckman et al1l2001j) and 
NGC 625 ^Cannon et alJl2005Tl . Mid-infrared images of 
the Galactic cen ter show a large-scale bipolar ou tflow in 
our own galaxy (Bland-Hawthorn & Cohen 2003). Anal- 
ysis of a larger sample of positive- velocity wings will be 
necessary to fully explore the properties of an O VI out- 
flow. 

9. SUMMARY 

We have studied FUSE and fiST/STIS spectra to 
investigate the ionization, kinematics, elemental abun- 
dances, and physical conditions in the HVCs toward 
HE 0226-4110 and PG 0953+414. We summarize our 
study in the following key points. 

1. HVCs toward HE 0226 4110. HVC absorp- 
tion in the range 80-230 km s -1 is detected in the 
HE 0226-4110 direction in lines of H I, C II, C III, 
C IV, O VI, Si II, Si III, and Si IV. Four com- 
ponents are clearly seen in this velocity range in 
the high-resolution STIS observations of C II, Si II, 
and Si III, at 99, 148, 175, and 193 kms" 1 . Each 
of these components is accompanied by a broader, 
highly ionized (C IV or O VI) component, indicat- 
ing that warm and hot phases exist in each compo- 
nent. The warm high- velocity clouds are all found 
to have similar densities, pressures, physical sizes, 
and ionization levels; their metallicities ([Z/H] be- 
tween —0.9 and —0.4) indicate the clouds may 
be outlying fragments of the Magellanic Stream, 
whose main body lies 11° away on the sky at sim- 
ilar velocities to the HE 0226-4110 absorbers. 

2. HVCs toward PG 0953+414. HVCs exist 
at both negative and positive velocities toward 
PG 0953+414. The negative-velocity HVC at 
-146kms _1 is seen in C II, Si II, Si III, and Si IV, 
has [Z/H]=— 0.6 ± 0.2, and contains no detectable 
C IV or O VI. At positive velocities, a narrow 
discrete HVC component with [Z/H] =-0.8 ± 0.2 
is seen in C II at 124 km s -1 . C III and Si III 
show much broader absorption components cen- 
tered at the same velocity. The highly ionized 
species Si IV and C IV show a weak component 
near to 150 km s -1 , offset from the lower ionization 
species by ~ 25 km s -1 . The O VI profile contains a 
smooth positive- velocity wing extending from 75- 
200 km s -1 . We show that part of this O VI absorp- 
tion could be produced in the same absorber that 
contains C IV and Si IV, but the gas in the range 
75-140 km s -1 requires an alternative explanation. 

3. CLOUDY Photoionization Modeling. The 

Galactic ionizing radiation field from hot stars 
dominates the EGB ionizing radiation inside the 
log < i> < gi2 = 4 contour, which occurs at approxi- 
mately 40 kpc for directions |6| < 30°, 100 kpc for 
30 < \b\ < 60°, and 180 kpc for directions |6| > 60°. 
Photoionization is a good explanation for the C II, 
C III, Si II, and Si III in the HVCs. Neither the 
EGB or the MW field can account for C IV and 
O VI in the HVCs, confirming that they exist in a 
separate, hotter gas phase. 
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4. Highly Ionized HVC Components. After 
studying the ionic column density ratio and kine- 
matic information of thirteen highly ionized HVC 
components reported here and in the literature, 
we believe the highly ionized HVC components 
arise at the interfaces between the cool/warm 
HVCs and a surrounding medium. Conductive 
interfaces can reproduce N(C IV)/iV(0 VI) and 
N(N V)/N(0 VI) in 11/12 highly ionized HVCs, 
though there are discrepancies remaining to be re- 
solved with the N(Si IV) prediction. HVC high-ion 
column density ratios are not signifi cantly d iffer- 
ent fro m the Galactic Halo averages of lZsargo et al.l 
( 2003), suggesting that similar ionization processes 
are at work in both environments. The interface 
scenario implies the existence of a hot, extended 
corona to the Milky Way extending out 50kpc to 
the Magellanic Clouds. Such a corona can also pro- 
vide confinement to the warm HVCs; our modelling 
shows the HE 0226-4110 HVCs would be close to 
pressure equilibrium in a fully ionized 10 6 K corona 
with riH = 4 - 9 x 10~ 5 cm~ 3 at 50kpc. 



5. Positive- Velocity Wings. The extended O VI 
absorption wing seen between 75 and 140 km s -1 
toward PC 0953+414 may trace a hot Galactic out- 
flow, because of the asymmetric profile and the lack 
of accompanying wing absorption in less-ionized 
species. The lack of C IV and Si IV in this velocity 
range can be used to constrain the temperature in 
the wing gas to be logT > 5.31, assuming CIE and 
solar elemental abundance ratios. 
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TABLE 1 
Sight Line Properties 
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z 


I 
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V 


-Fl030 a 










(°) 


(°) 


(mag.) 


(flux units) 


(kms -1 ) 
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Seyl 


0.495 


253.94 


-65.78 


15.2 


2.7 


-14.3 


PG 0953+414 


Seyl 


0.234 


179.79 


+51.71 


14.5 


5.1 


-0.3 



Note. — Information retreived from the NED database I http://nedwww.ipac.caltech.edu |, unless noted below. 

a Measured flux at 1030 A; 1 flux unit = 10~ 14 erg cm -2 s _1 A -1 . 

b Correction from heliocentric to LSR velocity frame. Av^gji = «lS-R — v Helio- 
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TABLE 2 
Summary of Observations 






Instrument 


Target 


Aperture 


Dataset ID 


Observation Date 


(ks) 


_&ST/STIS a 


HE 0226-4110 
PG 0953+414 


0.2" x 0.06" 
0.2" x 0.2" 


O6E107020.30 

O6E108010,20,30 

O6E109020,30 

O6E110010,20,30 

O6E111010,20 

O6E111030,40 

04X001010,20 
O4X002010 


2002 Dec 26 
2002 Dec 26 
2002 Dec 27 
2002 Dec 29 

2002 Dec 31 

2003 Jan 01 

1998 Dec 11 
1998 Dec 04 


6.0 
8.2 
6.0 
8.2 
5.1 
6.0 

Total=39.6 
13.7 
10.7 

Total=24.4 


FUSE b 


HE 0226-4110 
PG 0953+414 


30" x 30" 
30" x 30" 


P2071301 (v2.1.6) 
P1019101 (v2.1.6) 
P1019102 (v2.4.0) 
P1019103 (v2.4.0) 
P1019104 (v2.4.0) 
D0270101 (v2.4.0) 
D0270102 (v2.4.0) 
D0270103 (v2.4.0) 

P1012202 (v2.4.0) 
P1012201 (v2.1.6) 


2000 Dec 12 

2001 Oct 03 

2002 Nov 15 
2002 Nov 16 

2002 Nov 17 

2003 Sep 01 
2003 Sep 03 
2003 Oct 31 

1999 Dec 30 

2000 May 04 


11.0 
50.3 
14.5 
18.9 
18.1 
23.9 
39.9 
16.7 

Total=193.3 
38.9 
36.5 



Total=75.4 

a All STIS observations were taken with the E140M grating, providing 6.7kms -1 resolution (FWHM). The unbinned pixel size is between 
3.0 and 3.2 km s — ^varying with wavelength). 

b The FUSE resolution is 20kms -1 (FWHM) in the LiF channels, and M 25kms -1 in the SiC channels. Unbinned pixels are 2.0kms _1 
wide. We list the CALFUSE pipeline version used in parentheses after the Dataset ID. The difference between v2.1.6 and v2.4.0 is negligible. 




Fig. 1. — Maps in Galactic coordinates (7, b) of the HVC velocity field in the regions around HE 0226-4110 and PG 0953+414, based on 
the H I 21-cm data of Hulsbosch & Wakkcr ( 1988) and Morras et al. ( 2000). Colors show v^gji, and contours are at brightness temperature 
levels of 0.05, 0.5, and 1 K, corresponding to column densities of about 2, 20, and 40xl0 18 cm -2 . Individual cloud cores are labeled in the 
PG 0953+414 panel. 



TABLE 3 

Measurements of HVC Absorption toward HE 0226-4110 



Line 


S/N a 




HV1 


80 to 125 km s 


-l 




HV2, 125 to 160 k 


m s 1 




HV3 


, 160 to 185 


mi s 1 




HV 


4, 185 to i 








6 C 


W x d 


log N a c 


v b 





1/1/ d 


\ogN a c 




b c 


W x d 


\ogN a a 


v h 


b c 




H I A920.963 


10 


92±5 


25±2 


136+07+08 


1 P. OQ+0.02+0.04 
10.ZM_ 02-0 .04 


147±6 


14+2 


88+05+10 


1 r- 1 /-+0.06+0.18 

16.16+ J u J , 7 I ( ,. 3 i 


172+5 


12+1 


85+03+17 


1 P 9O+0.08+0.30 
10-zU_o ,n_o.30 


201+6 


19+3 


107+07: 


H I A919.351 


9 






f 


f 


146+7 


15+2 


79+05+06 


1 « 01 +0.05+0.05 
iD ' z± -0.07-0.06 


173+6 


12+1 


72+04+11 


i_ 9^+0.07+0.07 
10. zo_ .07-0.08 


203+6 


19+2 


110+07: 


H I A918.129 


8 


90±7 


27+2 


97+12+07 


1D.OU_ og _o .03 


147+7 


15+2 


64±07±09 


_i_n n74-ii nx 
iD - iO -0.08-0.08 


174+5 


11+1 


70+04+11 


1 P OCTU.uOT'J.UI 

10. oo_ .08-0.10 


199+6 


17+2 


95+07: 


H I A917.181 


8 


90±14 


29±3 


48±23±07 


i n i 9-Lfi on 
1D - uz -0. 17-0.07 


145+6 


15±2 


70+08+07 


_i_n nfi-Mi da 
^•"^-0. 07-0. 04 


173+5 


10+1 


68+05+11 


IP A A T^u.uo-f-u.uy 
1D.44_ .05-0 .11 


197+7 


15+3 


62+12: 


C II A1334.532 


9 


104±7 


16±3 


52±06±04 


lO.OZ_ 07-0 .03 
1 cn+0-05+0-06 
10. OU_ .06-0.06 


145+3 


13+1 


91+05+03 


io.yu_ 09 _o 01 


175+9 


10+3 


56+04+03 


- ~ „^4-0 ni-i-0 n r i 
1 _>.Dv>_q 05-0 .07 


198+5 


14+3 
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C IIA1036.337 


22 
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18±4 


33±04±07 


146+6 


15+2 


52+03+07 


10 77+0. 03+0. 06 
' ' -0.03-0.07 






f 


f 






f 


C III A977.020 
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105±8 


19±6 


83±06±09 


10. zo_ 02-0.05 


146±6 


12±6 
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- -. -.-.-1-0 (144-0 7(1 
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197+6 


15+6 
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C IVA1550.770 


9 






<45 


<13.35 


150+7 


14+5 


•37+1 n+i 1 

• ft _l L V ) _l L _L 


1 ^ qfj+0.14+0.05 
IO. OD_ .20-0.05 


167+7 


8+2 


25+09+08 


1 "i 9n+013+0.04 
lO.ZU_ .18-0.05 


213+6 


13+4 


35+12: 


NI A1200.710 


6 






<43 


<13.88 






<21 


<13.58 






<26 


<13.66 






<34 


N II A1083.994 


6 






<55 


<13.68 






<42 


<13.57 






<33 


< 13.46 






<40 


N V A1238.821 


10 






<44 


<13.32 






<21 


<13.00 






<20 


<12.97 




•g 


<47 


O I A1302.169 


12 






<23 


<13.48 






<29 


<13.58 






<16 


<13.33 




r-t- 

19*2 


<18 


O VIA1031.926 


22 


90+3 


33±2 


21+04+03 


n i n i n n"7 
-| q 9rTU.lUtU.Ur 

±o - z ' J -0.11-0.07 


145±6 


15+1 


26+03+03 


lO.OD_ .06-0 .09 


173+6 


11+1 


23+03+06 


-10 ni +u. Uo+U. UD 
IO. Ot_ .06-0.07 


201+7 


22+04: 


O VIA1037.617 


24 






<20 


<13.50 






i 


f 






f 


f 




•pr 


f 


Si IIA1260.422 


9 


102+7 


18±3 


47±08±04 


19 co+0.10+0.02 
iz.oo_ 12 -0 .01 


145±3 


1 9+9 
izztz 


RO-Lnfi-Lfl/l 

DZztUDltU4 


1 9 S1 +0-07+0.02 
iz.oi_ 08-0 .02 


174+8 


11+2 


42+05+03 


19 ko+0.07+0.04 
iz.oo_ .09-0.05 


199+6 


1*3 
18^5 


51+07: 


Si IIA1193.290 


6 


101±9 


21±6 


27+12+04 


19 7n+°- 20 + - u:i 

—0 37—0 03 


145±4 


12+3 


34+10+03 


1 9 7O+0.15+0.02 
— 23—0 01 






<42 


<12.76 


205+5 


25+12: 


Si IIA1190.416 


6 


105+6 


15±4 


32+13+03 


10 n c+0:i9+0.'03 
io.uo_ Q 35_o .03 


138+7 


10+3 


22+11+09 


19 oq+0.'18+0.'04 
iz.o»_ .29-0.04 


176+6 


8+3 


29+09+09 


19 qq+0.13+0.02 
iz.»»_o .20-0.03 






<60 


Si IIA1526.710 


10 






<17 


<12.88 






<42 


<13.27 






<14 


<12.80 


206+5 


18^4 


24+08: 


Si IIIA1206.500 


6 


88+7 


9+4 


25+08+08 


1 9 99+O.I6+O.O2 
^'"-0 .26-0 .02 


146+5 


14±2 


79+07+06 


19 7C+0.07+0.02 
iz. iu_ .08-0.02 


173+8 


10+1 


63+06+07 


1 9 C7+0-08+0.04 
iz - D '-0.09-0.03 


202+5 


13+3 


100+08: 


Si IVA1393.755 


10 






<34 


<12.57 






<30 


<12.52 






<21 


<12.36 


202+3 


10+5 


37+07: 


Si IVA1402.770 


12 






<21 


<12.67 






<23 


<12.72 






<37 


<12.92 


200+5 


18+2 


38+09: 


P II A1152.818 


15 






<34 


<13.09 






<25 


<12.95 






<23 


<12.92 






<16 


SII A1259.519 


10 






<34 


<14.17 






<32 


<14.14 






f 


f 






f 


S II A1253.811 


11 






<26 


<14.23 






<25 


<14.22 






<20 


<14.12 






<21 


S III A1012.495 


20 






<20 


<13.71 






<12 


<13.48 






<15 


<13.58 






<17 


Arl A1048.220 


26 






<16 


<12.81 






<12 


<12.68 






<18 


<12.86 






<14 


Fe II Al 144.938 


16 






<17 


<13.14 


145+2 


8+3 


11+4+5 


10 m +0.10+0.14 
lO-Ul_o .12-0.21 






<10 


<12.88 






<13 


FeIlA1608.451 


6 






<91 


<13.81 






<45 


<13.50 






<41 


<13.46 






<75 


FeUI A1122.524 


17 






<17 


<13.27 






f 


f 






f 


f 






f 



-Signal-to- 



,(«)'«/ Svt T a(")d" (kms" 1 ). 
: ^2/J+(v - v)^T a (v)dv/ /"+ T („)d«. (kms- 1 ). 



next to the lin 



Equivalent width measured between quoted velocity limits (mA), with statistical and systematic errors. 
C Apparent column density between quoted velocity limits, with statistical and systematic errors, or. for non-detections, 3cr limit . 
f 

No measurement or upper limit possible, because of blending with other lines. We list here the contaminated line, contaminated < 
HV3-4, C II* A1037.018 at v es 0; C IV A154S.195, HV1-4, N I A1200.223 at ^0; N I A1200.223, HV1-4, N I A1200.710 at u%0: O VIA1037 
Fe III A1122.524, HV2-4, IGM LyC at z = 0.20701. 



sity using 3cr equivalent width li: 



nd lii 



■ COG. N a 



the colun 

ncnts, blending line and velocity: H I A919.351, HV1, O Ia919.658 at V RiO; C II A1036.337, 
17, HV2-4, H 2 A103S.156 at u^O; S II A1259.519, HV3-4, Si II A1260.422 at u%0; 
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TABLE 4 

Measurements of HVC Absorption toward PG 0953+414 



Line 


S/N a 




HVl, - 


165 to -125 kins" 1 




HV2, 80 to 160 km s" 








v h 


b c 


Ws d 


iu g 1 v a 


v b 


b c 


"A 


log_/_ c 


tj T \oon orq 
rl 1 AyzU.yoo 


y 


1 rr n 1 n 

— 10o±o 


1 n 1 A 


89±07±10 


16 nq+ ' 02 + ' 05 
—0.03 — 0.07 


114±0 


23±3 


125±09±08 


1R 9R + 0.04+0.01 
lu.zu_ Q g_ 02 


H I A919.351 


9 


-149±6 


16±4 


74+07+- 1 n 


1fi ,0+0.04+0.07 
iu.lo_ Q 3 _o .07 






f 




H I A918.129 


9 


-138±8 


16±2 


77+08+-1 n 


1fi 9/.+0.03+0.07 
±Q.Z__ 03 _o .08 


109±7 


1 _+-8 


7Q^1 9_|_nn 


1 R on+O-OS + O- 02 
lu.zu_Q.io_o. 02 


H I A917.181 


9 






f 


f 


118±10 


Z _lEO 


£4+1 1 +-Oq 

o _^i izcuy 


1 R on+0.08+0.03 
lU._z_o 09— 0.02 


C IIA1334.532 


12 


-147±6 


12±4 


Z.A rir 1 no 


1 q c-7+0.05+0.02 
lo.O <_ .05-0.02 


116±6 


1 K 1 7 

lo± / 


At\ I on I no 

4y±uy±uo 


1 3 _q+0.07+0.01 
10. -0_q 08— 0.01 


C IIA1036.337 


12 






<38 


<13.53 






<46 


<13.59 


C III A977.020 


13 








f 


127±7 


40±5 


25±14±16 


i o 77+0.02+0.04 
lo. 1 1 _o 02—0 04 


C IV A1548.195 


9 






<38 


<12.97 


153±6 


23±4 


35±13±01 


12 QQ+0. 13+0.01 
—0.20—0.01 


C IV A1550.770 


8 






<47 


<13.36 






<38 


<13.27 


N I A1199.550 


7 






<45 


<13.42 






f 


f 


NI A1200.710 


7 






f 


f 






<86 


<14.19 


Nil A1083.994 


6 






<58 


<13.71 






<101 


<13.95 


N V A1238.821 


9 






<22 


<13.01 






<45 


<13.32 


O IA1302.169 


12 






<19 


<13.39 


113±10 


33±6 


23±08±02 


1 q cn + 0.12 + 0.03 

— U.l i — U. Uo 


O VI A1031.926 


14 






<39 


<13.49 


110±11 


37±5 


60±07±10 


1 S 7/1+0.05+0.06 
10. '*_o. 05-0. 06 


bi 11 AlzoU.4zz 


12 


— 147±7 


14±4 


44±04±02 


19 C7+0.05+0.01 

— U.Uo — U.U1 






<43 


<12.45 


o ■ TT \ 1 1 ah f_ r\ r\ 

bi ll All9o.290 


6 






<38 


<12.71 






<63 


<12.93 


Qi TT \1 1 on a i a 









<61 


<13.22 






<44 


/io n7 


Si IIA1526.710 


11 






<29 


<13.10 






<24 


<13.02 


Si IIIA1206.500 


9 


-144±2 


12±4 


55±07±08 


19 cc+0.03+0.06 
1Z.DO_ .02-0.07 


128±8 


42±4 


30±14±11 


19 on+o-os+o- 01 
iz.yu_ .05-0.01 


Si IVA1393.755 


11 


-144±7 


15±3 


20±07±03 


19 o Q +0.14+0.03 
iz.oy_ .i9_o.02 


150±6 


14±6 


18±08±02 


19 Q/+0.16+0.02 
iz.o__ .24-0.01 


Si IVA1402.770 


9 






<25 


<12.74 






<36 


<12.90 


P II A1152.818 


17 






<27 


<12.99 






<42 


<13.18 


S II A1259.519 


12 






<12 


<13.74 






<03 


<13.52 


S III A1012.495 


16 






<43 


<14.05 






f 


f 


Arl A1048.220 


18 






<26 


<13.02 






<26 


<13.02 


Fe IIA1144.938 


18 






<35 


<13.44 






<34 


<13.13 


Felll A1122.524 


18 






f 


f 






<47 


<13.72 



a Signal-to-noisc per resolution element in the continuum next to the line. 
b t. = Jj+ vr a (v)dv/ fv+ T a (v)dv (kms" 1 ). 

c b= t/2/.+(« - v) 2 T a (v)dv/ T a (v)dv. (kms" 1 ). 

^Equivalent width measured between quoted velocity limits (mA), with statistical and systematic errors. 

c Apparcnt column density between quoted velocity limits, with statistical and systematic errors, or, for non-dctcctions, 3cr limit set on the column 
density using 3- equivalent width limit and linear COG. N a is incm~ 2 . 

No measurement or upper limit possible, because of blending. We list here the contaminated line, contaminated components, blending line and 
velocity: H IA919.351, HV2, O IA919.658 at „Ri0; H IA917.181, HVl, edge of detector; C III A977.020, HVl, H 2 A976.552 at uf_0; N I A1199.550, 
HV2, N I A1200.223 at v ssO; N I A1200.223, HVl, N I A1199.550 at v a 0; N I A1200.223, HV2, N I A1200.710 at v f_ 0; N I A1200.710, HVl, 
N I A1200.223 at v ss 0; O VI A1037.617, HVl, H 2 A1037.146 at v f_ 0; O VIA1037.617, HV2, H 2 A1038.156 at v a 0; S III A1012.495, HV2, H 2 
A1012.822 at v «0; Fc III A1122.524, HVl, Fe IIA1121.975 at v «0. 

g Mcasurcments of C III, Si III, S III, Fe III, and O VI in HV2 are made over the extended velocity range 80 to 200 km s 1 ; measurements of 
C IV, N V, and Si IV are made over 125 to 200kms _1 . 



TABLE 5 



Summary of HVC Column Density Measurements 


Sight Line 


Comp. 


HI 


CII 


cm 


CIV 


Si II 


Si III 


Si IV 


O VI 


HE 0226-4110 
PG 0953+414 


HVl 
HV2 
HV3 
HV4 
HVl 
HV2 


16.29±0.05 
16.21±0.10 
16.74±0.20 a 
16.34±0.10 
16.15±0.15 
16.26±0.15 


IT CO+0.05 

lo.Oo_ 07 
13 qn +0 - 34 

io.yu_ 09 

,o fi o+0.06 
lo.uo_ 07 
1 S no+0.25 

io.yz_ o 5 

,o ^7+0.05 
lO.O(_ 05 

13.48_ Qg 


i_s._;8_ 05 
13.59l u; l° 
13.92_ ;°°3 
io.yu_ 05 

1 Q 77+0.25 
J-J' ' ' -0.04 


<13.35 

13.361°;^ 

13.20_ ;i* 

13 Sl+ - 17 
lo.ol_0.28 

<12.97 
12-99±8:_„ 


12 69+ ' 10 
lZ.OO_o.lo 

12.81±°;g| 
12 7n+ - 07 

1Z. (U_ 07 
19 71 +0.08 
L -0.08 
19 c 7 +0.05 
— 0.05 

<12.45 


1 9 99+O.I6 

iz.zz_ .2 6 

12 7fi+ - 07 
lz./D_ 08 

12fi7 +0.09 

LZ -°< -0.09 

12.89t°;g 7 

12 W+ a07 
iz.oo_o 07 

12.90±S:» 


<12.57 
<12.52 
<12.36 

12.85l°;« 

12-39±8.« 

12.34lg;« 


13 29+ ' 13 
io.zo_ 12 

13.36+0-1? 

,0 o 4 +0.08 
lo.o^_o 09 

13.29lg;?? 
<13.49 

,0 74 +0.08 
lo. ( __ 8 



NOTE. — All measurements in this table are the logarithm of the ionic column density in cm - . Errors are ±ler; upper/lower limits are 3er. 
For cases where absorption column densities are measured in more than one line (particularly Si II), we considered the dispersion of individual 
measurements and their errors to determine the final value listed. 

a Component saturated , even in weakest Lyman lines; column density estimated using AOD correction method for saturated lines 
ISavase fe Sembachll99ll) . 



Multi-phase HVCs 
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TABLE 6 

Flux of Ionizing Radiation near the Galaxy 





\j 


1V1UUL1 


l_slo lallLL 


lncrd) ,nm^ 
iu S ^<912 


lo g ^>912 




(°) 




(kpc) 


(<3> in cm -2 s _l ) 


($ in cm~ 2 s _1 ) 


HE 0226-4110 


-66 


MW 


10 


5.0 


7.0 








50 


4.5 


6.5 








100 


3.8 


5.8 






EGB 




4.0 


1.8 


PG 0953+414 


+52 


MW 


10 


5.5 


7.5 








50 


4.7 


6.6 








100 


3.8 


6.2 






EGB 




4.0 


4.8 



a MW = escaping Milky Way radiation (see §6.2); EGB = extragalactic background radiation (Haardt & Madau 1996). 
b Total ionizing flux impinging on the clouds at hydrogen-ionizing (hard; A < 912A) and non-hydrogen-ionizing (soft; 2460 A 
> A > 912 A) wavelengths. We assume F v = A-kJ v when computing $ for the EGB (see §6.1). 



TABLE 7 

Results of HVC Photoionization Modeling 



Sight Line 


Comp. 


V 


Model 


Dist. 


log[/ a 


[Z/H] b 


logn H c 


Z d 


P/k c 


/(H II) f 


7(Ha) s 






(km s" 1 ) 




(kpc) 






(nji in cm" 3 ) 


(pc) 


(cm" 3 K) 




(mR) 


HE 0226-4110 


HVC1 


99 


MW+EGB 


10 


-3.75 


-0.6 


-1.7 


11 


470 


0.98 


2.2 










50 


-3.70 


-0.6 


-2.1 


39 


170 


0.98 


0.8 










100 


-3.80 


-0.6 


-2.4 


56 


79 


0.97 


0.3 








EGB 




-3.80 


-0.5 


-2.6 


86 


51 


0.97 


0.2 




HVC2 


118 


MW+EGB 


10 


-3.55 


-0.4 


-1.9 


23 


290 


0.98 


1.9 










50 


-3.55 


-0.4 


-2.3 


57 


110 


0.98 


0.7 










100 


-3.55 


-0.3 


-2.7 


150 


47 


0.98 


0.3 








EGB 




-3.55 


-0.3 


-2.9 


230 


30 


0.98 


0.2 




HVC3 


175 


MW+EGB 


10 


-3.75 


-0.9 


-1.7 


37 


550 


0.98 


6.3 










50 


-3.70 


-0.9 


-2.1 


120 


190 


0.98 


2.3 










100 


-3.70 


-0.9 


-2.5 


290 


78 


0.98 


0.9 








EGB 




-3.70 


-0.9 


-2.7 


440 


50 


0.98 


0.5 




HVC4 


193 


MW+EGB 


10 


-3.50 


-0.6 


-1.9 


38 


240 


0.98 


2.6 










50 


-3.50 


-0.6 


-2.3 


98 


95 


0.98 


1.0 










100 


-3.50 


-0.5 


-2.7 


240 


39 


0.98 


0.1 








EGB 




-3.45 


-0.4 


-3.0 


500 


25 


0.99 


0.2 


PG 0953+414 


HVC1 


-150 


MW+EGB 


10 


-3.70 


-0.6 


-1.2 


3 


1180 


0.98 


5.2 










50 


-3.65 


-0.6 


-2.0 


21 


200 


0.98 


0.9 










100 


-3.66 


-0.6 


-2.5 


81 


59 


0.98 


0.3 








EGB 




-3.60 


-0.5 


-2.8 


170 


36 


0.98 


0.1 




HVC2 


-125 


MW+EGB 


10 


-3.15 


-0.9 


-1.8 


56 


460 


1.00 


5.0 










50 


-3.15 


-0.8 


-2.5 


310 


86 


1.00 


0.9 










100 


-3.11 


-0.8 


-3.1 


1290 


24 


1.00 


0.2 








EGB 




-3.20 


-0.7 


-3.2 


1270 


18 


0.99 


0.1 



Note. — Our photoionization models arc conducted by finding the values of the ionization parameter U and metallicity [Z /H] that best reproduce 
the observed ionic column densities of H I, C II, C III, Si II, and Si III, for a given input radiation field at a given distance. 



a Bcst fit ionization parameter, ±0.2 dcx (95% c.l.). 

b Best fit metallicity, ±0.2dex (95% c.l.), or ±0.3dcx for HE 0226-4110 HVC3 and HVC4, where saturation may affect the C III line. 

c Total hydrogen density: nu — n-y/£/cm~ 3 , where n-y is the density of ionizing photons in the radiation field. 

^Cloud depth: / — -/Vh/^h, where Nu is the total hydrogen column density returned by the model. 

G HVC pressure, P/k — ^ . n^T, summing over all elements and ionization states i and their associated electrons. 

^Hydrogen ionization fraction. 

g Predicted intensity of Ha emission: Z(Ha)=2.75(T/ 10 4 K)°' 924 n^ Rayleighs (R), where 1 R = 10 6 /4-rr photons cm -2 s" 1 sr~\ 
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TABLE 8 

High Ion Column Density Ratios in HVCs a 



Sight Line 


V 

(kms" 1 ) 


JV(H I) 
JV(0 VI) 


JV(Si IV) 
JV(0 VI) 


JV(C IV) 
JV(0 VI) 


AT(N V) 
AT(0 VI) 


JV(C IV) 
AT(Si IV) 


HE 0226-4110 
PG 0953+414 


+ 99 
+148 
+175 
+193 
-150 
+124 


nuu_ 400 

71 0+ 240 
' lu -230 

2500+ 1500 
zouu_ 1100 

1 1 no+ 400 

iiuu_ 300 
>600 

oo +150 
OJU -120 


<0.12 
<0.08 
<0.05 
0.36 ± 0.16 
>0.11 
0.04 + 0.02 


<0.72 
1.0 + 0.5 
0.7 + 0.3 
1.0 + 0.5 

0.18 + 0.07 


<0.80 
<0.26 
<0.23 
<0.65 

<0.19 


>10 
>10 

2.9 + 1.8 
<2.6 
4.5 + 2.5 


PG 1259+593 b 
Mrk 279 c 
PKS 2155-304 d 

Mrk 509 d 

PG 1116+215° 


-110 
-140 
-140 
-170 
-240 
-300 
+100 
+184 


(1.7±0.3)xlO b 
(6.8 ± 1.6) xlO 5 
370+ifo 
47+66 

-19 

< 2700 

< 2600 

1350 ± 40 

7000+ 3000 
(uuu_ 2000 


0.10 + 0.02 
0.22 ± 0.07 
0.09 ± 0.02 
0.05 ± 0.02 
< 0.02 
0.4 + 0.2 
0.18 + 0.06 
0.12 ± 0.01 


0.35 ± 0.05 
0.89 + 0.29 
0.51 + 0.06 
1.1 ± 0.2 
0.44 ± 0.07 
1.7 + 0.4 
1.12 + 0.36 
0.55 + 0.07 


< 0.07 
0.19 + 0.07 

< 0.08 

< 0.08 

< 0.07 

< 0.09 

< 0.27 

< 0.05 


3.4 + 0.4 

4.0 + 1.6 
5.9 + 1.1 

1.1 ± 0.2 
> 25 

4.2 + 2.0 

6.3 + 1.3 

4.5 + 0.8 


HVC average (detections only)' 
Halo average 5 






0.17 + 0.13 
0.20 + 0.13 


0.80 + 0.42 
0.60 + 0.47 


0.19 + 0.07 
0.12 + 0.07 


4.1 ± 1.6 
3.5 + 1.1 



a Refer to Tables 3 and 4 for the velocity ranges and log N a measurements for each component. Errors given are ±lcr; upper and lower 
limits are 3o\ 
b From lFm et alJ 120041) . 

c New measurements of N(C IV) and N(Si IV) using archival spectra of Mrk 279. We measure N(C IV)=13.68iq ^ioo^ and 
N(Si IV)=13.0ll° [Uto ° \ in the range v LS r=-210 to -115kms _1 , where high-velocity O VI is detected. 
dFrom lTollins et alJ ll2p0^ 
e From [Ganguly et al l 120051) . 

f Mean value of detections , with quoted error representing sample standard deviation. Upper/lower limits are ignored. 
s From Zsareo et al. (2003). Errors represent sample standard deviation. 



TABLE 9 

Origin Models for Highly Ionized HVC Components 



Origin 


Ionization 
Source 


Explains 3 
N(Si IV) , 
JV(0 VI) • 


Explains 3 
JV(C IV) 7 
N(0 VI) • 


Explains 3 
JV(N V) ? 
N(0 VI) ' 


Ref. 


Galactic Fountain 


Recombination of hot gas 


0/8 (4/5) 


1/12 (1/1) 


0/1 (11/12) 


(1) 


HVC interface 


Conductive Heating 


2/8 (4/5) 


11/12 (1/1) 


1/1 (11/12) 


,2) 




Turbulent Mixing 


7/8 (2/5) 


6/12 (0/1) 


1/1 ( 6/12) 


(3) 




Shock Heating 


6/8 (5/5) 


11/12 (1/1) 


0/1 (12/12) 


(4) 



References. — (1): lEdgar fc Chevalier! l ti986T) : (2): IBorkowski eTaT] Il990h : (3): ISlavin et alHl993T) : (4): IDopita k, Sutherland! IT996T) . 

a In these columns we assess what fraction of highly ionized HVCs ratios (from this paper and the literature) can be reproduced (to within 
their ±1<t errors) by the given model assuming sola r abundances. Th e numbers in parentheses represent the fraction of upper/lower limits 
that are consistent with the model predictions. See lFox et all 120041) for a tabulation of ionic ratio predictions in different models. 
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Fig. 2. — Absorption line profiles along the line of sight to HE 0226-4110. Flux is plotted against LSR velocity with the continuum 
and zero levels shown as dotted lines. The upper left panel shows the 21-cm emission line profile. Absorption in the HVCs is shaded, and 
blends are annotated at the expected line center of the offending feature. For O VI A1037.617 and C II A1036.337, the continuum includes 
a model of the H2 lines. (N) indicates night-only data is displayed to reduce geocoronal airglow emission. 
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Fig. 3. — H I Lyman series absorption line profiles along the line of sight to HE 0226-4110. Dotted lines display the position of the 
continuum and the zero level. Geocoronal H I emission centered near Okrns - 1 is seen in the stronger lines and can be ignored. The red 
line shows a Voigt component model of the H I absorption, using the component velocities and column densities derived using the AOD 
method. 
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Fig. 4. — Comparison of apparent column density as a function of velocity for various species seen in the HVCs toward HE 0226—4110. 
The units on the y-axes are ions cm -2 (kms -1 ) -1 . Error bars showing the typical uncertainty on each data point are shown in the corner 
of each panel, and small red arrows indicate where the high-velocity absorption may be saturated. The color of the shading indicates the 
velocity range covered by each component. 
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Fig. 5.— Same as Figure 2, except for PG 0953+414. 
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Fig. 6. — Same as Figure 3, except for PG 0953+414. The dotted red line shows the Voigt profile model without a component at 
-150 km s -1 ; this model clearly does not match the data, showing that H I exists and can be measured in the — 150kms _1 HVC. 
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Fig. 7. — Same as Figure 4, except for PG 0953+414. The color scheme reflects the shape of the absorption profile: narrow (blue), broad 
(green), broad and offset (red), and wing-like (yellow). The thick solid line in the O VI panel shows the predicted contribution from a 
conductive interface with a velocity and width defined by the C IV and Si IV profiles, and with a typical C IV/O VI ratio. 
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Fig. 8. — Radiation fields incident upon the HVCs i n our CLOUDY c a lculat ions. The blue line shows the estimated extragalactic 
background flux, where F v = AmJ v and J„ is taken from Haardt & Madau (1996). The green line shows our new estimate of the flux of 
escaping Galactic radiation, at distances of 0, 10, 50, and 100 kpc along the PC 0953+414 sight line (b = +52°). Note the substantial drop 
at 54 eV in the Galactic spectrum, caused by the He II edge in hot stars. Other ionization edges are marked at the top of the plot. 
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Fig. 9. — Contour plot showing the flux of Galactic radiation as a function of location relative to the Galactic disk in our new disk-halo 

model. The blue contours show the lines of constant hydrogen ionizing flux ($012 = Jq 12 (F\/hu)d\), and the red contours shows the 

(-2460, 



lines of constant non-hydrogen ionizing flux (<J > >gi2 = /gj*2 " (F\ / hv)dX) . The numbers next to the contours are the logarithm of <J><gi2 in 
units of cm -2 s — x . The Galactic ionizing field dominates the extragalactic ionizing background inside the line defined by log<I> = 4. The 
dot-dashed lines show the sight lines of interest in this paper (PG 0953+414 above the plane and HE 0226-4110 below it). This map is 
used to scale the normalization of the Galactic disk ionizing spectrum shown in Figure 8 for use in our photoionization modeling. 
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Fig. 10. — The behavior of various column density ratios of carbon and silicon as a function of ionization parameter, for HVCs immersed 
in three radiation fields: the EGB, the MW (any distance), and MW+EGB at 50 kpc toward HE 0226-4110. The logarithmic differences 
in the ionic ratios between the pure MW and pure EGB fields are shown in the panels labelled MW-EGB. Under the assumption of pure 
photoionization, these curves can be used to convert an observed ionic ratio to an implied ionization parameter, for a given radiation field. 
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Fig. 11. — Example CLOUDY run showing the fits to observed column densities in HVC2 toward PG 0953+414, for two radiation fields: 
EGB (left) and MW+EGB at lOkpc (right). For both radiation fields a good fit can be found to the C II, C III, Si II, and Si III data, but 
not the C IV, Si IV, or O VI. The best-fit log U and [Z /H] and their errors are derived using a chi-squared minimization process, shown 
in the grayscale on the bottom panel. The contours show the regions of 70, 95, and 99% confidence. Results from models of this kind to 
all HVC components toward HE 0226-4110 and PG 0953+414 arc given in Table 7. 
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Fig, 12, — High -ion column density ratios measured in HVCs (red points, purple upper limits), the Galactic halo average (green point; 
IZsargo et al. 2003), and the predictions of various models (gray boxes). CIs = conductive interfaces; CIE = collisional ionization equilibrium; 
TMLs = turbulent mixing layers; RC = radiative cooling; SI = shock ionization. The bl ue arrows show t he effect of correcting the model 
predictions with solar abundance ratios onto a Magellanic Stream set of abundances. See lFox et al.l 120041) for a detailed explanation of the 
model parameters. 



